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I ABSTRACT 

, For a sample of long Gamma-Ray Bursts (GRBs) with known rcdshift, wc study 

CIh' the distribution of the evolutionary tracks on the rest-frame luminosity-peak energy 

Q I LisQ-Sp diagram. We are interested in exploring the extension of the 'Yonetoku' cor- 

^ . relation to any phase of the prompt light curve, and in verifying how the high-signal 

' prompt duration time, Tj, in the rest frame correlates with the residuals of such 

Ci . correlation (Firmani et al. 2006). For our purpose, we analyse separately two sam- 
ples of time-resolved spectra corresponding to 32 GRBs with peak fluxes Fp > 1.8 

^ ! phot cm-2 s~i from the Swift-BAT detector, and 7 bright GRBs from the CGRO- 

^ ■ BATSE detector previously processed by Kaneko et al. (2006). After constructing the 

OO \ Liso~Ep diagram, we discuss the relevance of selection effects, finding that they could 

. affect significantly the correlation. However, we find that these effects are much less 

•/^ ' significant in the LisoTj-E'^ diagram, where the intrinsic scatter reduces significantly. 

" I , We apply further corrections in order to reduce the intrinsic scatter even more. For 

7—i • the sub-samples of GRBs (7 from Swift and 5 from CGRO) with measured jet break 

T— I I time, tj, we analyse the effects of correcting Ljso by jet coUimation. We find that (i) 

OO , the scatter around the correlation is reduced, and (ii) this scatter is dominated by 

' the internal scatter of the individual evolutionary tracks. These results suggest that 

^ \ the time integrated 'Amati' and 'Ghirlanda' correlations are consequences of the time 

• i resolved features, not of selection effects, and therefore call for a physical origin. We 

' finally remark the relevance of looking inside the nature of the evolutionary tracks. 

H : 

^ Key words: gamma rays: bursts - gamma rays: observations 



1 INTRODUCTION 

The phenomenological study of long Gamma-Ray Bursts 
(GRBs), on one hand, brings the possibility to use these ex- 
treme cosmic events as tracers of many astronomical and 
cosmological processes, and on the other, leads to a compre- 
hensive description of the underlying GRB physics (see for 
recent reviews Meszaros 2006; Ghirlanda, Ghisellini & Fir- 
mani 2006; Zhang 2007). The finding of correlations among 
spectral and energetic properties of the prompt emission has 
been the main avenue of progress in these directions. 

Amati et al. (2002) discovered a correlation between 
the isotropic-equivalent energy radiated during the whole 
prompt phase, Eiao, and the peak energy of the rest-frame 
time-integrated ufv spectrum, < E'^ >. The scatter around 
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the Eiso-< Ep > correlation suggests the existence of some 
hidden parameter(s). Ghirlanda, Ghisellini & Lazzati (2004) 
have identified such parameter with the jet coUimation an- 
gle, 6j, and Liang & Zhang (2005) with the rest-frame 
break-time in the afterglow light curve, tj ; 6j is tightly asso- 
ciated to tj. Both parameters are time integrated quantities 
of the outfiow in the fireball model. A correlation analogous 
to the Amati' one but for the peak isotropic-equivalent 
luminosity, I/iso,p, has been first reported by Yonetoku et 
al. (2004). Later on, for a sample of 19 GRBs, Firmani 
et al. (2006) have found a strong correlation among I/iso.p, 
< E'p > and a rest -frame prompt emission high-signal GRB 
time duration, Tj . Rossi et al. (2008) have also found a cor- 
relation among I/iso.p, < E'p >, and T^, but with diff'erent 
slopes and scatter. More recently, CoUazzi & Schaefer (2008) 
confirmed the 'Firmani' correlation, though the scatter they 
obtained for a larger sample is higher than in Firmani et al. 
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(2006). The homogeneity and the quahty of the data seems 
to be crucial for obtaining reliable statistical results; in this 
sense, the publicly available large database from the Swift is 
by now an invaluable source of data. 

The GRB light curves are composed by a sequential su- 
perposition of many pulses, possibly originated by individual 
shock episodes. The spectral and temporal characteristics 
of these pulses are key ingredients for understanding the 
prompt emission mechanisms of GRBs (e.g., Ryde & Pet- 
rosian 2002). Although there is not a simple pattern to de- 
scribe the great variety of the light curve morphologies, some 
relations have been found in GRBs with light curves charac- 
terized by a single pulse or a few well-separated pulses. For 
example, Borgonovo & Ryde (2001) have found a correla- 
tion between the (bolometric or spectral peak) energy flux, 
F, and observed Ep during the decaying phase of individual, 
well-defined pulses: F oc -Ep'', with a mean 7 « 2. Notice 
that Ep in this case is inferred from each time-resolved spec- 
trum. Further, Liang, Dai & Xu (2004) have tested whether 
or not the relation F oc i?p^ within a burst holds in a sub- 
sample of 2408 time-resolved spectra for 91 bright CGRO- 
BATSE GRBs presented by Preece et al. (2000; all kind 
of pulses are included in this sample). They reported that 
for 75% of the bursts, the Spearman correlation coefficient 
is larger than 0.5. Nevertheless, even for the Borgonovo & 
Ryde sample of well-defined pulses, this correlation has a 
considerable scatter, which is probably produced by the su- 
perpositions of pulses. 

If the individual pulses within a burst follow on average 
the correlation F oc Ep'^, then a natural question is whether 
there exists or not a universal relationship among all GRBs 
in the rest frame. 

Both the 'Yonetoku' and the Liso,p-< E'p >~T'j-- correla- 
tions (Firmani et al. 2006) were established for GRBs with 
known z taken from heterogeneous samples, with incom- 
plete available data, and by using hybrid quantities: the fiux 
used to calculate Liso.p is at the peak of the light curve, 
while both < Ep > and the bolometric correction had to be 
inferred from the time-integrated spectrum. The primary 
data publicly available from the Swift satellite (Gehrels et 
al. 2004) offers now the possibility to study the mentioned 
correlations in the rest-frame and with temporal resolution. 
Unfortunately, the narrow spectral range of the Swift BAT 
detector limits and makes difficult the spectral analysis. 

By means of time-resolved spectral analysis for a sam- 
ple of 32 Swift and 7 bright CGRO long GRBs, we wiU 
show here that, when multiplying Liso by T'f, the evolution- 
ary tracks of different GRBs form a well defined correlations, 
with a small scatter. In §2 we develop a heuristic framework 
to introduce such a connection. The GRB samples used and 
the spectral analysis are presented in §3. The construction of 
correlations with the prompt time-resolved data, the anal- 
ysis of the selection effects, and a discussion of the results 
are presented in §4.1. In §4.2, we present our attempt to 
further reduce the scatter around the correlations by intro- 
ducing the jet coUimation correction. We show that most of 
the remaining scatter in the overall correlations is due to 
the individual scatter of each evolutionary track. In §5 we 
present our conclusions. 




90 
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Figure 1. Maximum bolometric flux of a given event versus 
arctan(7) for 208 CGRO GRBs with 5 or more time-resolved 
spectra (from the Kaneko et al. 2006 sample). The quantity 7 is 
the slope of the best linear fit to the F~Ep evolutionary tracks. 
The arctan function is used in order the error bars (Sc) looked 
symmetric. The vertical dashed line indicates arctan 2 = 63.4 de- 
grees. Thin red error bars show the events with 7 = 2 within 3cr, 
while thick black error bars indicate the few events incompatible 
with 7 = 2 at the 3a level. There is not any systematical trend 
of ^max with 7, however, the uncertainty and the spread of 7 
around 2 increases as -Fmax decreases. 

2 MOTIVATIONS 

Our purpose in this paper is to study, in the rest frame, 
evolutionary features of the prompt 7— ray spectra and their 
connection with other GRB properties. Before discussing the 
data, we develop an heuristic reasoning aimed to show a 
potential connection between the phenomenological time- 
integrated and time-resolved correlations mentioned in the 
Introduction. 

We start from the fact that the prompt evolution of a 
given burst can be characterised, at least above some flux 
and/or some peak energy limit, by the dependence F oc Ep'' , 
where F and Ep are defined at any given time. The power 
index 7 as well as the scatter of the evolutionary correlation 
could change from burst to burst. Previous analysis showed 
that the mean value of 7 is around 2 with some dispersion 
(Borgonovo & Ryde 2001; Liang et al. 2004). In the rest 
frame, the burst evolutionary track F oc Ep'' is equivalent 
to the correlation 

Li,o = KE'p\ (1) 

In the logarithmic plane E'p-Liso the evolutionary track de- 
scribed by Eq. ((T} identifies a straight strip with a slope 7, 
a zero point, and a specific sector on the strip. 

Let us first assume 7 = 2. Note that K in this case is 
independent of the Doppler-boosting factor D = l/[r(l — 
/3/i)], where F, /3c, and fi are the Lorentz factor, the ve- 
locity, and the cosine of the viewing angle, respectively. By 
integrating in time, the evolutionary track reduces to 

Ei,o = KT' <E'p^>, (2) 

where T' is a rest-frame time scale of the energy emission. 
Now, if observations give a correlation between Eiso and 
< E'p >^ among different GRBs, (the 'Amati' relation), then 
KT' — k becomes approximately a constant independent 
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of each GRB. The conclusion is that scaling K oc T'^^ , 
the scatter of the correlation given by Eq. ((T)) for a sample 
of GRBs is less than the scatter around the I/iso-< E'p > 
correlation for the same sample. The fact that 7 = 2 for 
all the GRBs simplifies the previous reasoning: since the 
evolutionary track given by Eq. and the 'Amati' relation 
show the same slope, then only the zero point of the strip 
is involved in the scatter around the 'Amati' relation. The 
parallel strips superimpose and thus reduce their occupation 
region on the plane. 

If 7 7^ 2, the situation becomes more intriguing. Now 
the evolutionary tracks are not aligned with the 'Amati' re- 
lation and so it is more difficult for them to reduce their 
occupation region on the E'p-Liso diagram. Nevertheless, if 
not only the zero points of the strips, but even the occu- 
pation sectors of the evolutionary tracks within the strips 
conspire in order to agree with the 'Amati' relation, then 
a reduction of their occupation region on the plane (reduc- 
tion of the scatter) may be reached. This inference allows to 
hope that T' plays a role of a hidden parameter related to 
the scatter of the correlation Eq. ([T|. 

This heuristic reasoning has actually inspired the 
present work. Definitively, it is of great interest to verify with 
high quality prompt time-resolved data whether or not the 
'Firmani'-like time-resolved relation E'^ vs LisoT' is obeyed, 
and to study the scatter around it. If this relation is obeyed, 
then a corollary is that it establishes a link between the 
'Amati' _Eiso — and the 'Yonetoku' Liso^p — < iJp > 

relations and, even more importantly, a robust physical in- 
gredient has to exist behind such relation in the sense that at 
each instant of the prompt the GRB ts aware of the duration 
of the entire process. 

Before presenting results from our spectral analysis on 
Stmft and CGRO GRBs with known redshift, we explore the 
large set of time-resolved spectra from the CGRO data anal- 
ysed by Kaneko et al. (2006). Since for the great majority 
of the bursts in this sample the redshift is unknown, we are 
able to infer only the values of the slope 7 of the individual 
correlations F cx E^'' . The slope should remain the same for 
the corresponding individual rest-frame correlations Liso vs 
E',. 

The Kaneko et al. (2006) sample contains 350 bright 
GRBs and 7427 time-resolved spectra. Each time-resolved 
spectrum was fitted to a Band (Band et al. 1993) function. 
We consider only time-resolved spectra with Band model 
parameters a > —2; /3 < — 2 and 50<i5p<1800 keV (note 
that Ep is in the observer frame) with an uncertainty better 
than 50% and a time step shorter than 4 s. We assume a 
bolometric correction given by the Band energy distribution. 
Finally, we take into account only GRBs with at least 5 
time-resolved spectra. Such conditions reduce the sample 
to 208 GRBs. 

For this sample, only 16% of the GRB tracks show a 
disagreement with 7 = 2, while the other 84% are compati- 
bles with 7 = 2 at the 3(t level. A similar result, but for less 
events than here, has been reported by Liang et al. (2004; 
see their figure 2 for the F — Ep tracks corresponding to 
some bursts). We plot in Fig. [1] the maximum bolometric 
flux, Fmax, of a given event versus arctan(7); arctan is used 
in order the error bars (at 3cr) looked symmetric. It is clearly 
seen that only a small fraction of the events (thick black er- 
ror bars) have values of 7 different from 2 within 3(j (vertical 



dashed line shows arctan 2 = 63.4 degrees); the majority of 
them (thin red error bars) show 7 compatible with 2 within 
the 3(t level. In Fig. [T] it is also seen that the correlation 
slope 7 is not biased systematically by Fmax, though its un- 
certainty and the spread around 2 increse as -Fmax decreases. 
The result shown above opens the question of identifying the 
physical mechanisms that determine the value of 7 as well 
as the intrinsic scatter of each evolutionary track. The fact 
that a large fraction of GRBs shows 7 ~ 2, (i) supports 
the consistency of our research concerning a universal cor- 
relation, (ii) provides a natural GRB identification criterion 
to eventually optimise such correlation, and (in) allows to 
study the behaviour of those GRBs with 7 7^ 2. 



3 SAMPLE SELECTION AND SPECTRAL 
ANALYSIS 

The time-resolved spectral analysis requires data with a 
high enough signal-to-noise ratio in order to ensure a rea- 
sonable determination of the spectral parameters. We select 
here 32 long GRBs from the whole Swift-BAT sample of 
GRBs with known redshift (until February 2008) and with 
a peak fiux, _Fp, greater than 1.8 phot cm~^ s~^. Because 
of the Swift-BAT narrow effective spectral range (15-150 
keV), a careful statistical handling of the correlated uncer- 
tainties in the spectral parameters is mandatory (Cabrera 
et al. 2007). Even so, a significant fraction of the observed 
GRBs and of the time-resolved spectra within a given GRB, 
had to be discarded because the peak in their ufi, spectra 
lies out of the BAT limit or because the signal is too low. 
At the end, we have 32 Swift usable bursts with 207 time- 
resolved spectra. 

The photon model adopted to fit the Swift time- 
resolved spectra is the cut-off power law (GPL), which has 
three parameters. The fits are carried out with the heasoft 
package XPSEC 0. The time-resolved spectra selected for 
the analysis are chosen to be shorter than the light curve 
variation time scales, but large enough as to ensure accept- 
able confidence levels (CLs) for the photon index and Ep. 
We reject the prompt time-resolved spectra where Ep can 
not be determined or the uncertainty on Ep and/or on the 
photon index is too large due to the scarce number of counts. 
With these constraints, we obtained 207 usable time slices 
for our 32 Swift GRBs. Table 1 reports the basic information 
concerning this sample: the name, the redshift, the observer- 
frame Tf and jet break time tj, and finally the rest-frame 
time-integrated E'p and total Eiso- Tf is the time spanned in 
the observer frame by the brighest 50% of the total counts 
above background (Reichart et al. 2001) for the light curve 
on the observed energy range 17-100 keVB A more detailed 
description of the temporal binning and selection criteria 
will be presented elsewhere. Liso is calculated from the GPL 
spectral parameters within the energy range 1-10000 keV 
at rest. The correlated errors of the spectral parameters are 
adequately propagated in order to obtain the corresponding 

^ http: / /heasarc. gsfc.nasa.gov/docs/softwarc/lheasoft/download. html 

We measure Tf as in Reichart et al. (2001), but instead of 
using the duration of the brightest bins in the light curve that 
enclose / = 45% of the total counts, we have used / = 50%. 
None of the results presented here changes by such redefinition. 
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Figure 2. Bp vs Liso for the time— resolved spectra corresponding to 32 Swift GRBs (68% CL ellipses) and 7 bright CGRO GRBs (Icr 
error bars). The solid lines connect the data-points of each individual GRB (evolutionary tracks); red colour is for the Swift data and 
green colour for the CGRO data. The best linear fit to the Swift data is shown with the solid right line, while the ±lcr intrinsic scatter 
is indicated with the parallel dot-dashed lines. The upper solid curve corresponds to a limiting observed {z = 0) flux of 10~^erg cm~^ 
s~^ and E-p ~ 100 keV. The ticks along the curve indicate different ^'s: 1, 2,. ..,9, from left to right. 



correlated errors (CL ellipses) of Liso and E'p (see for details 
Cabrera et al. 2007). 

With the aim to compare it with a completely diflter- 
ent sample, along with the 32 GRBs observed by Swift, we 
have included in our considerations 7 bright CG-RO-BATSE 
GRBs with known redshifts reported in Table 2. Now Tf is 
measured on hght curves in the 50 — 300 keV observed energy 
range. For these 7 GRBs we have 287 useful time-resolved 
spectra available from Kaneko et al. (2006). In the case of 
the CGRO sample, the spectral information tends to be of 
much better quality and the time-resolved spectra can be 
fitted with the more general four-parameter Band model. 

The Swtft-BAT and CGRO-BATSE samples studied 
here are different in many aspects. We remark two of them. 
The first one is due to the spectral model, CPL for Swift 
and Band for CGRO. Based on Cabrera et al. (2007), a 
rough estimate of this difference may be obtained adopting 
for the Swift spectral fitting a Band model with /3 frozen to 
—2.3. The result shows an increase in Liso by a factor 1.6 for 
a given E'p. The second difference is the light-curve energy 
range in which Tf is estimated. For Swift we use 17—100 keV, 
while for CGRO, we use 50 — 300 keV. Taking into account 
the inverse dependency of on the energy at the power 0.4 
given by Reichart et al. (2001), the factor 3 between Swift 
and CGRO energy ranges leads to a Tf for the Swift sample 
roughly 3°'^ = 1.55 larger than T} for the CGRO sample. 



Curiously enough, both effects leave the product Liso x Tf 
roughly invariant. In spite of this coincidence, we have de- 
cided to handle each one of the samples separately, and only 
afterwards we will take care to compare the results of our 
analysis that result invariant with respect to the systematic 
differences mentioned above. 



4 CORRELATION ANALYSIS 

4.1 Prompt 7— ray emission correlations 

Figure [5] shows in the Liso-E'p diagram the time-resolved 
spectral data for the Swift sample (ellipses) and for the 
CGRO sample (error bars). The ellipses correspond to the 
68% CL error regions calculated taking into account the 
error covariance matrix (see Cabrera et al. 2007), while the 
orthogonal error bars correspond to the standard deviations. 
No correction for the differences in the spectral model has 
been applied (see §2). The scatter in the plot is rather large. 
The red lines show the evolutionary tracks for the Swift sam- 
ple. The brightest (highest S/N ratios) bursts show 7 ~ 2 
while the fainter bursts show values between 1 and 4. It is 
not possible to obtain more continuous and detailed tracks 
because of relevant sections of the light curves that have 
time-resolved spectra with Ep lying outside the BAT spec- 
tral sensitivity. The green lines show the evolutionary tracks 
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name 


z 


(sec) 


(days) 


Lg < Ep' > 
(keV) 


(10^" erg) 


N 


050315 


1.949 


15.00(0.50] 








1 


050318 


1.440 


4.00(0.10] 




2.45(0.04] 


2.07(0.11] 


2 


050401 


2.900 


6 20f0 101 


1 50 fO 501 


3 45 fO 231 


3 51 fO 1 41 


3 


050505 


4.270 


11.80(0.10] 




3.54(0.16] 


3.13(0.12] 


2 


050525 


0.606 


3.10(0.02] 


0.30(0.10] 


2.33(0.01] 


2.34(0.01] 


19 


050603 


2.821 


2.10(0.10] 




3.68(0.30] 


3.55(0.22] 


3 


050820A 


2.612 


12.00(0.40] 


15.00(8.00] 


3.12(0.09] 


3.99(0.03] 


2 


050922C 


2.198 


1.44(0.02] 




3.36(0.22] 


2.63(0.14] 


6 


051 1 oq 


2.346 


6 90 fn 201 




2 97f0 1 61 


2 80f0 1 91 


\ 


051111 


1.549 


1 1 60f0 101 




3 36 fO 341 


2 80f0 221 


1 


0fi020fi 


4.048 


2 70f0 1 01 




3 29f0 061 


2 59 fO 061 


5 


ofi02i n 


3.910 


35 00f3 001 


2 20f0 901 


S 45f0 1 41 


3 52 fO 101 


4 


060418 


1.490 


16 20 fO 301 








3 


060526 


3.221 


12.30(0.30] 


2.80(0.30] 


2.02 (0.09] 


2.41(0.05] 


1 


06061 4 


0.125 


28 70 fO 1 01 


1 38 fO 041 


1 74 fO 361 


1 40f0 1 71 


g 


060904B 


0.703 


10.00(0.30] 




2.36(0.13] 


1.44(0.13] 


5 


060906 


3.685 


11 50 fO 201 




2 99f0 091 


3 IfifO 141 


X 


060927 


5.600 


4.10(0.20] 




3.50(0.06] 


2.92(0.06] 


4 


061007 


1.262 


18.40(0.08] 








5 


061121 


1.314 


6 20f0 051 








3 


07031 8 


0.836 


12 00 fO 701 








9 


070508 


0.820 


6.10(0.05] 








20 


070521 


0.553 


10.30(0.20] 




2.66(0.10] 


2.11(0.07] 


22 


070810A 


2.170 


3.30(0.10] 




2.58(0.07] 


2.14(0.15] 


5 


071003 


1.100 


12.80(0.50] 








9 


071010B 


0.947 


4.42(0.05] 




2.25(0.04] 


2.29(0.09] 


14 


071117 


1.331 


1.50(0.10] 




2.99(0.21] 


2.21(0.12] 


7 


080319B 


0.937 


21.11(0.07] 


9.00(2.00] 


3.12(0.01] 


4.12(0.01] 


2 


080319C 


1.950 


5.16(0.08] 




3.15(0.15] 


2.63(0.11] 


5 


080411 


1.030 


6.55(0.02] 








21 


080413A 


2.433 


6.37(0.03] 








9 


080413B 


1.100 


1.15(0.01] 




2.59(0.12] 


2.19(0.10] 


7 



Table 1. Sample of 5Mi/t-BAT long GRBs with known z and useful time-resolved spectra. The prompt observer emission time duration 
Tf was calculated by us (see text), while the break time tj was taken from the literature compilation by Ghirlanda et al. (2007); for 
GRB080319B, was taken from Kann, Schulzc & Updike (2008). The time— integrated rest— frame peak energy and isotropic-equivalent 
energy were taken from Cabrera et al. (2007) or calculated by us using the GPL model. N is the number of useful time— resolved spectra. 



name 


z 




h 


Lg < Ep' > 


LgBiso 


N 






(sec) 


(days) 


(keV) 


(lO-M erg) 




970828 


0.957 


11.00(2.00] 


2.20(0.40] 


3.06(0.09] 


3.47(0.05] 


35 


980703 


0.966 


18.80(0.30] 


3.40(0.50] 


2.99(0.09] 


2.84(0.05] 


2 


990123 


1.600 


19.50(0.10] 


2.00(0.50] 


3.72(0.03] 


4.38(0.05] 


84 


990506 


1.307 


15.30(0.20] 








80 


990510 


1.619 


5.50(0.10] 


1.60(0.20] 


3.05(0.04] 


3.25(0.05] 


20 


991216 


1.020 


4.30(0.05] 


1.20(0.40] 


3.12(0.09] 


3.83(0.05] 


54 


000131 


4.500 


2.00(0.50] 








12 



Table 2. Sample of CGi?0-BATSE long GRBs with known z and useful time— resolved spectra calculated by Kaneko et al. (2006). 
In this case, the Band model was used for both the time-resolved and time-integrated quantities. See the caption of Table 1 for more 
details. 



for the CGRO sample. Four events show 7 = 2 at la, two 
(990123 and 990510) at 2a, and one (990506) at 3a. We have 
shown in §2 that even for the much larger sample of CGRO 
GRBs without redshift determination (from Kaneko et al. 
2006), indeed 7 ~ 2 for most of the events. 

The bottom right part in the Liso-E'p diagram is physi- 
cally empty; here, the selection effects do not apply. There- 
fore, a real upper limits Liso for each E'p does exist. On the 



contrary, in the top left part of the diagram, selection ef- 
fects could be present. The upper continuous curve plotted 
in Fig.[5]corresponds to a limiting hypothetical observed flux 
of 10~*erg cm~^ s~^ between 15-150 keV and an observed 
peak energy ~ 100 keV. This curve gives a rough idea of the 
sensitivity limit of the Swift-'BAT instrument. Therefore, 
the present data do not allow to identify a specific boundary 
here because of the limiting fluxes characterizing the sam- 
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Figure 3. Panel (a): Liso/E'p^ vs T| for each burst from the 
same Swift time-resolved spectral data shown in Fig. [2] and for 
9 short Swift GRBs. The size represents the redshift while the 
colour represents the peak flux in phot cm~^ s~^. The dashed 
lines correspond to a —1 slope. Panel (b): Residuals 5 of the 
correlation Ljao-Ep shown in Fig. [2] and their corresponding un- 
certainties vs T'j. The best fit slope — 0.30±0.04 leads to a scatter 
reduction term proportional to j'j.i-25±0.20^ 



pies. Thus, taking care of eventual selection effects in this 
region of the diagram, a kind of 'Yonetoku' correlation may 
be extended to the time-resolved features and may include 
a conspicuous fraction of the prompt evolution. This corre- 
lation might be reflecting an intrinsic local physical process 
of the GRB emission mechanism. 

Along the constant fiux/i?p curve plotted in Fig. [2] we 
show with ticks the values that would have Ep and Liao at 
different redshifts (integer z from 1 to 9). This result shows 
that this correlation would not be useful to infer pseudo- 
redshifts for GRBs with non determined redshift. A similar 
conclusion will apply for the E'p-LiaoT'f correlation presented 
below. 

The thick straight solid line in Fig.[2]is the best linear fit 
to the logarithmic data performed by a minimum method 



taking into account the residual and its uncertainty on the 
orthogonal direction of the straight line. In Fig. [2] (as well 
as in the following related Figures), in order to avoid over- 
crowding, the best fit is plotted only for the Swift sample. 
The fit gives a line described by: 



^°s(i^)=«Log 



10'' erg s~ 



(3) 



The corresponding coefficients of the fits to both the Swift 
and CGRO samples are given in Tables 3 and 4, respectively 
(relation LE). Notice that h is not a parameter of the fit, but 
is the logarithm of Liso in the barycentre of the data-points. 
The standard errors were computed in the barycentre frame 
of the data-points in order to minimise any correlation be- 
tween them. The average correlation slope of both samples is 
0.30 ±0.05, and the difference in the intercepts at E'^ = 100 
keV between the CGRO and Swift best fits is 0.09 ± 0.11. 
We have also estimated the intrinsic scatter aintr of the cor- 
relation by adding it iteratively in quadrature to the error 
in the y axis (Logi?p) and by requiring that the reduced Xr 
be equal to 1. This method has been tested by Novak et 
al. (2006) and it gives values similar to those obtained by 
the fitting method presented in D'Agostini (2005). The in- 
trinsic scatters (standard deviations) around the E'^ vs Liso 
correlation for both samples are given in the same Tables 
3 and 4. The average value of the (Tintr is 0.21 ± 0.01. The 
dot-dashed lines in Fig. [2] show the intrinsic scatter around 
the Swift sample correlation. 

Taking into account the result by Firmani et al. (2006) , 
our next step is to explore a possible correlation of Liso and 
Ep with the high-signal emission rest frame time T'j. For 
calculating T'f, we have taken into account the corrections 
due to the cosmological time dilation and the narrowing of 
the light curve's temporal substructure at higher energies 
(Fenimore et al. 1995). Following Reichart et al. (2001), this 
last correction on Tf, calculated fixing the energy band on 
the rest frame, goes as (1 -I- zY with e ~ 0.4. Then apply- 
ing both corrections we obtain T'j = T/(l + z)^^'^'^ Actually 
the value of e is not well constrained. For a large BATSE 
database, Zhang et al. (2007) found that the standard devi- 
ation of the distribution of e (modeled as Gaussian) is 0.51 
with a median value of 0.39; the distribution is skewed to 
smaller values. 

For the Swift sample we find that Liso roughly tends to 
be oc {T'f)~^ with a large scatter, while E'^ does not show 
any significant correlation with Now, if we plot LisoZ-Bp"* 
from each time-resolved spectrum vs T'f, then the scatter 
reduces for m ~ 2 and a clear trend with slope —1 (upper 
dashed line) is observed (panel (a) of Fig. [3]). This result 
implies that the residuals in the Liso-E'p diagram (Fig. [2} 
correlates with T'j. 

It is interesting to show in panel (a) of Fig. [3] the po- 
sition of some short GRBs with known redshift. We have 
carried out time resolved spectral analysis for 9 Swift short 
GRBs. As seen in Fig. (3] these data are naturally segregated 
in the diagram, having significantly smaller Liso/E'p^ values 
than those of the long GRBs. 

With the aim of checking for possible selection effects 
behind the correlation presented in panel (a) of Fig. |3] we 
plot the data with the color code representing the observed 
peak photon flux (Fp) range, and with the cross size indi- 
cating the redshift range. Our first test has to do with the 
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relation a b c fintr 

LE 0.25±0.02 51.48 2.27±0.02 0.195±0.013 

LTE 0.39±0.02 52.17 2.27±0.01 0.137±0.010 

LTHE 0.49±0.03 49.48 2.11±0.02 0.093±0.020 

LTWE 0.66±0.05 49.42 2.14±0.02 0.105±0.020 



^tot 
intr 



0.084±0.007 
0.097±0.008 



Table 3. Best— fit slopes a and zero points c of the different logarithmic linear relations studied in this paper for the Swift sample. The 
parameter b defines the absisa of the barycentre (10*', in cgs units) of the given correlation while (Tintr is the intrinsic scatter. For the 
correlations which include the collimation angle, o-^°l^ is an estimate of the scatter contribution due to the internal dispersion of the 
evolutionary tracks. 



relation 



a 



Tintr 



LE 0.36±0.03 52.91 2.87±0.01 0.217±0.010 

LTE 0.42±0.02 53.72 2.88±0.01 0.174±0.008 

LTHE 0.50±0.02 51.23 2.94±0.01 0.104±0.007 

LTWE 0.53±0.02 50.67 2.93±0.01 0.131±0.008 



0.130±0.007 
0.130±0.007 



Table 4. Same as in Table 3 but for the CGRO sample. 



influence of _Fp on the Liso/E'p^ vs T[ diagram. We have 
divided the sample into three sub-samples: 1.8 < -Fp < 5 
(red), 5 < Fp < 10 (green), and > 10 (blue) (the units 
are phot cm~^ s~^). As seen in Fig. [S] the correlation is re- 
markably insensitive to the peak flux. In other words the 
sensitivity limit does not influence the correlation; it only 
limits the number of objects on it. Our second test refers to 
a possible bias with z. According to the same panel of Fig. 
[31 it is evident that the data plotted in the T'f-Li^o/E^^ dia- 
gram are not appreciably biased by z. A third test concerns 
a possible selection effect on the GRB duration. By plotting 
-Fp vs T'f for all the Swift GRBs with known z (even those 
with Fp < 1.8 phot cm s , which were not used in our 
analysis here), we have seen that the limit Fp > 1.8 is well 
above any biasing limit for Tf. 

Finally, while the upper right side of the diagram in 
panel (a) of Fig. [3] involves fluxes that are not biased by 
selection effects, the presence of short GRBs on the bottom 
left side of the same diagram gives a further evidence against 
any selection effect here. Unfortunately any further exhaus- 
tive analysis on the distribution of Fp is impossible due to 
the limited BAT spectral capability. We conclude that the 
correlation L\so/ E'^ vs r[ is reasonably free from selection 
effects. 

Panel (b) of Fig.Oshows for the Swift sample how much 
the residuals of the Liso-E'p diagram of Fig. [2] correlate with 
Ty. Here S are the orthogonal residuals of the Liso-E'p best fit 
(positive 5 correspond to high Liso), its standard deviation 
being calculated in the same direction. The best fit slope 
of the 5 vs T'j correlation is —0.30 ± 0.04. Combining this 
result with the Liso-Fp best fit slope we obtain that the 
correlation Fp vs Li^oTf^ reaches its minimum scatter for 
p — 1.28 ± 0.20. A similar analysis on the CGRO sample 
leads to p = 1.12 ± 0.25. Given the relevance of this result 
we have made use of other more sofisticated methods based 
on multilinear analysis obtaining in each case p = 1.25±0.20 
and p = 1.00 ± 0.20, respectively. This result implies that 
the emission time T'f of the events on the diagram of Fig. 
[2] increases along the orthogonal direction of the best fit 



straight line when one goes from the low Fp - high List, 
to the high Fp- low Liso- Therefore, as it has been found 
in Firmani et al. (2006) for a different sample, T'j reduces 
the scatter around the correlation Fp vs Fiso- Later on we 
will estimate the reduction of such scatter. This result is 
particularly intriguing because it reveals how instantaneous 
features such as Liso and Fp are actually regulated according 
to the overall duration of the burst. 
We fit the data with the line 

where we adopt p — 1.25 and p = 1 for the Swtft and CGRO 
samples, respectively. The best linear fit parameters to Eq. 
([4]) for the Swift and CGRO samples are given in Tables 
3 and 4, respectively (relation LTE). The average of both 
slopes is 0.40 ± 0.02, while the difference in the intercepts 
at Fp = 100 keV between the CGRO and Swtft best fits is 
0.15 ± 0.07. The average intrinsic scatter is (Jintr ~ 0.15 ± 
0.02. A remarkable decrease on the scatter from the Fp vs 
Fiso to the Fp vs LisoT'f'' correlation is evident. Figure [4] 
shows Fp vs LisoT'j'' for the same samples plotted in Fig. 
[21 The best fit line, as in Fig. [51 refers only to the Swift 
sample, and the corresponding scatter is represented by the 
dot-dashed lines. Our basic considerations will not change 
appreciably assuming p = 1 even for the Swift sample. In 
fact, the internal scatter changes from 0.137 to 0.142 with a 
standard deviation of 0.010. The p = 1 assumption will be 
taken later on just for economy. 

The previous discussion about selection effects on the 
diagrams of Fig. [3l can be translated now into the Fp vs 
LisoT'j-^ correlation (Fig. [4|. We conclude then that the lat- 
ter correlation is reasonably free of selection effects. This is 
rather evident if we imagine that the effect of is to shift 
the borders of the Fp vs Liso correlation into the body of the 
Fp vs LisoF/^ correlation. Then the lack of low luminosity 
events infiuences now the population of the Fp vs Li^oT'^^ 
correlation and not its borders. 
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Figure 4. Ep vs L\soT'j for the same time-resolved spectral data shown in Fig. [2] The symbol and line codes are as in Fig. [2] The 
best-fit line and the ilc intrinsic scatter (solid and dot-dashed lines, respectively) are shown only for the Swift sample. 



4.2 Scatter reduction by correcting luminosities 
due to collimation 

Following Ghirlanda et al. (2004) , a way to reduce the scat- 
ter around the E'^ vs LisoT/'' correlation could be by cor- 
recting the time-resolved isotropic luminosities by the jet 
collimation angle in order to estimate the intrinsic 7— ray 
temporal luminosities, L-y = (1 — cos9j)LiBo- The collima- 
tion semi-aperture angle 6j is calculated from the jet break 
time tj and the total emitted isotropic energy Eiso by two 
alternative models: the homogeneous ISM model (HM) and 
the wind medium model (WM). Unfortunately, reliable esti- 
mates of tj for our GRB sample are available only for some 
events. We use the tj values compiled from the literature by 
Ghirlanda et al. (20070. Our GRBs with known reduces 
to 7 events of Swift with 37 time resolved spectra, and 5 of 
CGRO with 195 time resolved spectra (see Tables 1 and 2). 
In the HM model, the jet opening angle is given by: 



^ Note that, in principle, should be achromatic, since it is due 
to a geometrical effect. However, we rarely have true achromatic 
breaks in the well sampled light curves of Swift bursts. This may 
be due to the fact that the optical and the X-ray emission are due 
to two different components (see e.g. Uhm &; Beloborodov 2007; 
Genet, Daigne & Mochkovitch 2007; Ghisellini et al. 2007). As 
discussed in Ghirlanda et al. (2007), it is likely that the optical 
emission is more often associated to the forward shock of the 
fireball running into the circumburst medium, and therefore more 
indicative of possible jet breaks. 



where r]^ is the radiative efficiency that we assume equal to 
0.2, and n is the medium density that we assume equal to 3 
cm"^. 

Figure[5]shows the collimation-corrected correlation E'^ 
vs. L-yT'f for 12 GRB with known tj (left correlation), to- 
gether with the same data presented in Fig. |4] (right correla- 
tion). The best fit line showed for the collimation-corrected 
correlation refers only to the Swift sample; the correspond- 
ing scatter is represented by the dot-dashed lines. Fitting 
the data with the line 

the best fit parameters of the E'^ vs L-^T'^ correlations (HM 
case) for Swift and CGRO samples are given in Tables 3 and 
4, respectively (relation LTHE). The slopes of both correla- 
tions are very close, the average being 0.50 ±0.03. Concern- 
ing the scatter, for both samples (Tintr has clearly reduced 
its value with respect to the one in the E'^ vs LisoT/ corre- 
lations; on average we find Uintr = 0.10 ± 0.01. 

We conclude that the HM collimation angle introduces 
a remarkable reduction on the scatter in the E'^ vs -LisoT/ 
correlation. However, such a reduction does not concern the 
internal scatter in each GRB evolutionary track. 

In order to estimate the contribution of the internal 
scatter from each evolutionary track to the overall correla- 
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Figure 5. The data and lines in the right part of the diagram are the same as those presented in Fig. [2] (the LisoT^--Ep diagram). The 
data in the left part correspond to the GRBs of both Swift and CGRO samples with known tj. For such GRB, has been used instead 
of Lisa (the HM case was used to calculate 9j). The symbol and line codes are the same as in Fig. |2] Note how the intrinsic scatter has 
been reduced by going from the L^^^T'^-E'^ diagram to the L-yT'^-E^ one. 



tion intrinsic scatter, we have performed the following exer- 
cise. For each sample (Swift or CGRO), every LisoT'j-E'p 
GRB track is shifted rigidly to best fit with the coUimation- 
corrected correlation E'p vs i-yTj presented above for the 
corresponding sample (see Fig. |6] for the Swift sample; here 
the evolutionary tracks of the CGRO GRBs are included 
only for graphic display). We are now able to calculate 
the new intrinsic scatter of the corresponding sample, cr'"',., 
around the respective correlation E'p vs L-yT'f. The values 
of o-,'°tj (relation LTHE) are reported in Tables 3 and 4 for 
the Swift and CGRO samples, respectively. By comparing 
ffintr and o-fntr, it is rather evident that the scatter of the 
tight coUimation-corrected correlation is dominated by the 
internal scatter of the evolutionary tracks. A similar result 
is obtained if instead of the entire sample we estimate the 
a*°tj. taking into account the GRB with the known jet break 
time alone. In fact in this case the average internal scatter 
is o"intr = 0.08 ± 0.01. We conclude that any further reduc- 
tion of the coUimation-corrected correlation scatter could be 
reached by reducing the internal scatter of each evolutionary 
track, and that the latter should be possibly identifying the 
hidden parameters behind the stochastic properties of the 
evolutionary tracks of each GRB. 

In the WM model the jet opening angle is given by 



where now the medium density is supposed to be n (r) = 

5 X 10^^^*r~^ g cm~^ and we assume A, = 1. Tables 3 
and 4 present the corresponding best fit parameters for the 
Bp vs T'f correlation in the WM case (relation LTWE). 
In this case the best fit slopes are moderately close (95% 
CL). The average of both slopes is 0.60 ± 0.05. The scatters 
(Tintr of the WM coUimation-corrected correlations are also 
smaller than the ones of the not corrected correlations E'^ vs 
LisoT'f, but the scatter reduction is smaller than in the HM 
case. On average we find (Tintr = 0.12 ± 0.01. After perform- 
ing the same track shifting procedure described above, we 
find that also in this case the internal scatter of each GRB 
evolutionary track provides the dominant component of the 
scatter around the tight E'p vs L^tj correlation. 

Finally, we have explored whether the residuals of the 
LisoTf-E'p correlation correlate or not with several prompt 
light-curve parameters: the variability V (Reichart et al. 
2001); the "emission symmetry" Sf = Tb/Ti, where Ti and 
T2 are the duration times of the fiuence-halves, from 5 — 50% 
and 50 — 95% of the total counts, respectively (Borgonovo 

6 Bjornsson 2006); the Ep'-'^^ / Ep^^^ ratio, where Ep'^^^ and 
Ep^^'' are the peak energies of the integrated ufi, spectra for 
each of the two time intervals Ti and T2, respectively (Bor- 
gonovo & Bjornsson 2006). Our preliminary results show 
that the residuals are not correlated with any of these pa- 
rameters, i.e. none of them could be a potential reductor of 
the scatter around the LisoTf-Ep correlation. 
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Figure 6. The data and lines in tlie riglit part of the diagram are the same as presented in Fig. |4] (the Li^aT'^-E'^ diagram). The data 
on the left correspond to the rigid shift of each evolutionary track to best fit with the E'^ vs L~^T'j correlation for the Swift sub-sample 
of GRBs with known (left part of Fig.[5ll. The CGRO data are shown only for illustrative purposes; Note that the scatter of the E'^ vs 
LryT'r corrclatiou constructed by the track shifts is only slightly smaller than the one of the original coUimation-corrected correlation. 



5 CONCLUSIONS 

We have selected the high-signal time-resolved spectra from 
the available sample of Swift long GRBs with measured z, 
and analysed them with the aim to search for systematic 
features of the local 7— ray emission mechanism and their 
connection with known global GRB properties. Requiring a 
GRB peak flux Fp > 1.8 phot cm~^ s~^, a total of 207 time- 
resolved spectra corresponding to 32 GRBs (until February 
2008) were analysed. We have included also 287 spectra from 
7 bright CGRO GRBs with z known analysed previously in 
Kaneko et al. (2006). Since the two samples are affected in 
a different way by some systematic effects, we preferred to 
perform our correlation analysis separately for each sample 
and then to check whether the results are consistent or not 
between them. We have found that they are indeed con- 
sistent. Thus, for simplicity, in what follows we report the 
averages of the two samples for the best-flt parameters of 
the different correlations. 

The main results and conclusions from our study are as 
follows: 

• By plotting the time-resolved data-points in the loga- 
rithmic Liso~E'p diagram, a linear band with average slope 
0.30 ± 0.05 and intrinsic scatter aintr ~ 0.21 ± 0.01 appears. 
While the low Ep - high Liao region is free of selection ef- 
fects, the high E'p- low Liso region could be affected by the 
flux limits of the samples. 



• We found that the residuals in the Liso-E'p diagram cor- 
relate with T|. This result offers a strong evidence that the 
parameter T| reduces the scatter of the Ep vs Liso correla- 
tion. By analyzing the vs Liso/E'p'^ diagram (Fig. |3}, we 
have checked that selection effects are not responsible for 
such a trend. 

• In agreement with the previous point, we have intro- 
duced the logarithmic diagram E'p vs LiaoT'j^ and have found 
that the optimal value p = 1.1 ± 0.1 reduces the scatter to 
Cintr = 0.15 ± 0.02. The average slope of the correlation is 
0.40±0.02. Such correlation reveals three important aspects. 
First, its intrinsic scatter is smaller than the intrinsic scat- 
ter in the E'p vs Liao correlation. Second, it is reasonably 
free from selection effects, accordying to our our analysis in 
the T'j vs Liao/E'p^ diagram. Third, it represents a connec- 
tion between instantaneous features {E'p, Liao) and global 
features (T'f). At any moment, the instantaneous features of 
a GRB correlate with the entire prompt duration as if at 
each instant of the prompt the GRB would be aware of the 
duration of the entire process. 

• For the 12 GRBs out of our samples (7 from the Swift 
and 5 from the CGRO samples, respectively) for which the 
jet break time tj is known, we could further reduce the scat- 
ter of each sample by using the coUimation-corrected lumi- 
nosity L-y by estimating the collimation jet angle for the ho- 
mogeneous (HM) and wind medium (WM) cases. The lowest 
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intrinsic scatter has been obtained for the E'p vs L^T^ corre- 
lation in the HM case; the (average) slope of the correlation 
is 0.50 ± 0.03 and aintr = 0.10 ± 0.01. 

• We have estimated the contribution of the internal scat- 
ter of the evolutionary tracks to the scatter of the over- 
all correlations. For this, the Liao-Ep evolutionary track of 
each GRB has been shifted to best fit with the correspond- 
ing coUimation-corrected (HM and WM cases) correlations. 
Our results indicate that for both cases crintr ~ t^lntu this 
means that the total intrinsic scatter is mainly due to the 
internal scatter of the tracks. Thus, with the caveat that 
the statistics is still limited, we conclude that any further 
reduction of the scatter around the GRB empirical corre- 
lations may be attained by discovering the hidden variables 
behind the stochastic features of the individual Bp -Liso evo- 
lutionary tracks. 

We conclude that the long GRB individual evolutionary 
tracks populate a rather narrow strip in the Ep vs LiaoT/ 
diagram with a slope « 0.4, whatever the evolutionary track 
slope is. The jet coUimation correction further reduces the 
thickness of such strip and leads the slope close to 0.5. While 
selection effects probably are present in the E'p vs List, di- 
agram, they do not seem to weaken our conclusion. This 
implies the existence of a universal 7-ray emission mecha- 
nism for long GRBs where the instantaneous features are 
modulated by a global parameter, which we found here to 
be T^. We suggest that the interconnection between T'j, and 
the Liao-Ep evolutionary tracks is at the basis of the global 
'Amati' and 'Ghirlanda' relations. 
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